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Abstract. We consider a simple model of a cloud chamber consisting of a test particle (the 
a-particle) interacting with two other particles (the atoms of the vapour) subject to attractive 
potentials centered in ai,a 2 G M 3 . At time zero the a-particle is described by an outgoing 
spherical wave centered in the origin and the atoms are in their ground state. We show that, 
under suitable assumptions on the physical parameters of the system and up to second order in 
perturbation theory, the probability that both atoms are ionized is negligible unless a 2 lies on 
the line joining the origin with a±. The work is a fully time-dependent version of the original 
analysis proposed by Mott in 1929. 



1. Introduction 

The classical limit of quantum mechanics is a widely studied subject in mathematical physics 
and many detailed results on the asymptotic regime K — > are available (see e.g. [Rj and refer- 
ences therein). Nevertheless, in most cases, the limit h — > for the time-dependent Schrodinger 
equation relative to a given quantum system is studied only for suitable chosen, "almost clas- 
sical" states, namely WKB or coherent states. Roughly speaking, these results guarantee that 
if one choses an almost classical initial state then for h — > its propagation remains close to a 
classical propagation at any (not too long) later time. The problem arises when one considers 
a situation in which the quantum system at initial time is in a genuine quantum state, e.g. 
a superposition state, and nevertheless the system exhibits a classical behavior. Examples of 
this situation are the localization effect in chiral molecules or the suppression of the interfer- 
ence fringes for a heavy particle in a two-slit experiment. It is clear that the emergence of 
such classical behavior cannot be understood if one insists to consider the limit h — > for 
the isolated quantum system. It is worth mentioning that the problem has some relevance 
from the conceptual point of view. In fact it was already raised in the earliest debate on the 
foundation of Quantum Mechanics (see e.g. [BEJ). The accepted explanation in these cases is 
based on the consideration of the interaction of the quantum system with an environment, and 
in particular on the decoherence effect produced by the environment. It is important to stress 
that the decoherence effect must be proved in each situation, starting from specific models of 
system + environment and introducing suitable assumptions on the parameters of the model. 

l 
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Many results in this direction have been obtained in the physical literature (see e.g. the reviews 
[GJKKSZ] , [Horn] and references therein) but only few mathematical results are available. 
Here we want to focus on an old problem of a different kind, raised by Mott ((Mj) in 1929, 
concerning the explanation of the straight tracks left by an a-particle in a cloud chamber. Ac- 
cording to quantum mechanics ( [Gaj . [CGj ) the a-particle, isotropically emitted by a radioactive 
source, is initially described by a spherical wave function and then interacts with the atoms of 
the vapour, whose positions are randomly distributed in the chamber. The observed tracks are 
the macroscopic manifestation of the ionization of the atoms induced by the a-particle and "it 
is a little difficult to picture how it is that an outgoing spherical wave can produce a straight 
track; we think intuitively that it should ionise atoms at random throughout space" ([M]). 
The explanation proposed by Mott was based on a simple model describing the a-particle in 
interaction with only two atoms. Exploiting time-independent perturbation arguments, he con- 
cluded that the probability that both atoms are ionized is negligible unless the two atoms and 
the center of the spherical wave lie on the same line. In such rather indirect sense, this would 
explain why we see straight tracks in the experiments. We will remark on this aspect in section 
3. Notice that so far we have spoken of tracks due to the decay of an a-particle. In order to 
be able to speak of a track of an a-particle one must use the reduction postulate which states 
that the initial spherical wave becomes localized at the place where the ionization takes place. 
We mention that the problem is also discussed in [H] and later in [B] , and some further elabo- 
rations on the subject can be found in [BPTJ, [BrJ, [CCF'J, [CLj . |Haj . [St]. 
In this paper we reconsider the three-particle model of a cloud chamber. Under suitable as- 
sumptions on the parameters of the model, which will be specified later, we give a proof of 
Mott's result through a fully time- dependent analysis and up to second order in perturbation 
theory. The method of the proof is rather elementary and it basically relies on stationary and 
non stationary phase arguments for the estimate of the oscillatory integrals appearing in the 
perturbative expansion. The work extends to the three-dimensional case the result obtained in 
|DFTj for the simpler one- dimensional case, where the spherical wave reduces to the coherent 
superposition of two wave packets with opposite mean average momentum. 



2. DESCRIPTION OF THE MODEL 

In this section we describe the model, i.e. the Hamiltonian, the initial state, the assumptions 
on the physical parameters, and introduce some notation. 

Let us introduce the Hamiltonian. We consider a three-particle non relativistic, spinless quan- 
tum system in dimension three, made of a particle with mass M (the a-particle) and two 
other particles with mass m which play the role of model-atoms. More precisely we describe 
such atoms as particles subject to an attractive point interaction placed at fixed positions 
a\ , a-i G M 3 , with |ai| < \a%\ (as we shall see later, the result is still valid if the point interac- 
tion is replaced by the Coulomb potential). Moreover we assume that the interaction between 
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the a-particle and each atom is given by a smooth two-body potential V. We denote by R the 
position coordinate of the a-particle and by r 1; r 2 the position coordinates of the two atoms. 
The Hamiltonian of the system in L 2 (IR 9 ) is formally written as 

H = H + XHt (2.1) 

H = K + K 1 + K 2 (2.2) 

H x = V{r\R - n)) + V(j'\R - r 2 )) (2.3) 
where K denotes the free Hamiltonian for the a-particle 

h 2 

K = A R , (2.4) 

A > is a coupling constant and Kj, j = 1,2, is the Schrodinger operator in L 2 (M 3 ) with an 
attractive point interaction of strength — (47T7) -1 , 7 > 0, placed at aj. The operator Kj is by 
definition a non trivial self adjoint extension of the free Hamiltonian restricted on C^°(1R 3 \{a,j}) 
(for further details see e.g. [AGH-KHJ). The spectrum is 



n 

and the proper and generalized eigenfunctions are respectively given by 



a p (Kj) = {E } , E = - 7T —^, a c {K j ) = a ac {K j ) = [0,oo) (2.5) 



= ^ Cd'Hr - %)) , (°(x) = -±= ^ (2.6) 

h(r,k) = e^0 o (7- 1 (r-a,), 7 fc), ef>°(x,y) = - — - ^ (2-7) 

The parameter 7 has the physical meaning of a scattering length and it characterizes the 
effective range of the point interaction. From (12.61) it is also clear that 7 is a measure of the 
linear dimension of the ground state, i.e. of the atoms. 

The unperturbed Hamiltonian Hq is obviously selfadjoint and bounded from below in L 2 (R 9 ) 
and moreover the smoothness assumption on the interaction potential V (see theorems 1, 2) 
guarantees that the perturbed Hamiltonian H is also selfadjoint and bounded from below on 
the same domain of H . In particular this implies that the evolution problem associated with 
the Hamiltonian H is well posed. 
We choose the initial state in the product form 

%(R, ri ,r 2 ) =V(fl)Ci(rOC2(r a ) (2.8) 
where Q is defined in (12.61) and ifi(R) is a spherical wave defined as follows. Let us consider a 
gaussian wave packet localized in space around the origin, with standard deviation 7 and mean 
momentum P$ > along the direction u G S 2 . Integrating over the unit sphere S 2 , one obtains 
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= ~^~3/2 f(^ lR ) I diL e ^ > f&) = ( 2 - 9 ) 



^7 3/ " " " ' Jsp 

where M is a normalization factor and e > is the dimensionless parameter 



h 

s = — 2.10 

"o7 



By an elementary integration we get 



R 



Air A/" p 1 

^(i2) = ^^ r Me- 1 7- 1 |i2|), AT= — - 2 (2.11) 



4vr 7 / 4 - 



We remark that the characteristic length 7 appears in the definition of the Hamiltonian as well 
as in the initial state in such a way that the range of the interaction between the a-particle 
and the atoms, the linear dimension of the atoms and the localization in space of the spherical 
wave are all of order 7. As it will become clear further on, this is a crucial ingredient for the 
proof of our result. 

We also notice that the initial state f)2.8p corresponds to the situation considered by Mott, i.e. 
the a-particle emitted as a outgoing spherical wave and the atoms in their ground state. 
Let us describe the hypotheses on the physical parameters of the model. We assume 

£<1 (2.12) 



^L = 0(e), j = 1,2 (2.13) 

171 

3-0(«) (2.14) 

Condition (I2.12p means that the wavelength HPq 1 associated to the initial state of the a-particle 
is much smaller than the linear dimension of the atoms and the range of the interaction, which 
means that we are in a semi-classical regime for the a-particle. In C I2. 131) we assume that |ai|, 
1 02] are macroscopic distances with respect to the characteristic length 7 and in (I2.14p we 
require that the mass ratio is small. Finally we tacitly assume 

Ao^ = 0( £ ) (2.15) 

where vq = PqM -1 . Condition fl 2 . 1 5 [) is necessary in order to make reasonable the application 
of our perturbative techniques, even if it is not strictly required for the proof of our results. The 
above assumptions (12.121) . (12.131) . (I2.14p have some relevant physical implications. In particular 
from (I2.12p . (12.141) one sees that the binding energy of the atoms is small compared to the 
kinetic energy of the a-particle 
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i -> . °( £ ) ( 2 - 16 ) 

Furthermore the assumptions (12.121) and ( 12. 13ft imply two relations among the characteristic 
times of the system which will be relevant in the what follows. In particular we define the flight 
times to the atoms of the a-particle 

r, = ^, j = 1,2 (2.17) 

v 

the characteristic "period" of the atoms 

T a = 27T-- = 4vr— L (2.18) 
l-C'ol n 

and the transit time of the a-particle in the region where the atom are localized 

T t = - (2.19) 

Then one has 

?± = l* = 0( e ), 3 = 1,2 (2.20) 
tj v \aj\ 

£ = -L"»- = 0( l) (2.2!) 
T a 47r m P 7 

i.e. the transit time T t is small with respect to the flight times Tj but it is comparable with the 
characteristic period of the atoms T a . This means that the a-particle can "see" the internal 
structure of the atoms. 

Let us introduce some notation to streamline the presentation. 

&j = T%, J = l,2 (2.22) 
\ a j\ 

uj(y) = ^(l + y 2 ), yeR 3 (2.23) 

ht hTj h 2 t . . , 

a = W bj = W' ^/w^ ' J = 1 ' 2 (2 - 24) 

h ^,y) = iA^ fdxe-^(x,y)(: (x), CyeM. 3 (2.25) 



(2tt) 3 / 2 

g{Z,v) = V{Z)h{Z,y) (2.26) 

where F denotes the Fourier transform of F. Moreover we denote by || ■ ||^p>™, p>l,fi,mEff, 
the standard weighted Sobolev norm in R 3 , and by <S(M 3 ) the Schwartz space of real function 
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in M 3 . Finally Ck denotes a positive numerical constant, depending on k G N and, possibly, on 
the dimensionless parameters (12.241) . 

It is important to notice that, for any fixed t > r 2 and yj G IR 3 , the quantities a, bj, Cj are all of 
order one. In fact it is sufficient to notice that 



h | a.2 1 t , h \aj\ M ( h \ I a-2 1 t 



•• n b< = „ 1 , Cj = — — - — u(yj) (2.27) 

P o7 7 r 2 ' J P 7 7 ' ™ V/W 7 r 2 VW V ; 

and to use (1212)1 . (12331) . (l2TT4j) . 

We conclude this section noticing that our results (i.e. theorems 1, 2 in the next section) are still 
valid if the point interaction Hamiltonians are replaced by the Hamiltonians with a Coulomb 
potential — e 2 |r — a-,] -1 , as model Hamiltonians for the atoms. In fact, if we denote by r the 
Bohr's radius, we can write the ground state energy and the corresponding eigenstate as 
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Si(r) = ^Z (r»\r- aj )), S°(x) = ^ e -W (2.29) 



Moreover the generalized eigenfunctions are 

$,(r, fc) = <5>°(r \r - aj ),r k), (-A x - -j-J $°(x, y) = y) (2.30) 

It is now sufficient to replace 7, E , Q, <pj by r , So, Sj, $j and it is easily checked that the 
proofs work exactly in the same way. 



3. RESULT 



We are now in a position to formulate our result. We are interested in the computation of the 
probability that both atoms are ionized at time t > 0. An exact computation obviously requires 
the complete knowledge of the state \I/(£) of the system, which is not available. Following the 
original strategy of Mott we shall limit to consider the second order approximation ^2^) of the 
state ^(t) which, iterating twice Duhamel's formula, is given by 



# 2 (£) = e 
if 2 (t) = ^0 



iUH Q 



4 fdh e^ H ° 
ft Jo 



ti 



*0 



(3.1) 



(3.2) 
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Therefore we shall study the probability that both atoms are ionized up to second order in 
perturbation theory, i.e. 



V(t) = dRdk 1 dk 2 



dn dr 2 0i (ri, fci)0 2 (r 2 , k 2 )4/ 2 (R, n,r 2 , t) 



(3.3) 

Our main result is the characterization of the ionization probability V(t) for a fixed time t > t 2 
and it is summarized in theorems 1, 2 below. In theorem 1 we consider the case in which a 2 
is not aligned with a± and the origin and we show that the ionization probability decays faster 
than any power of e. 



Theorem 1. Let us fix t > r 2 , \ai\ < \a 2 \, d\ ■ d 2 < 1 and let us assume (I2.12p . (12.1 3p . f)2.14p . 
V G 5(M 3 ). Then for any k G N there exists > such that 

-2k 



l,k 



\a 2 



1 - -. — r + (1 - di-d 2 ) 



_2k-2 



(3.4) 



The constant in (13 .4p . which will be specified during the proof, remains finite and strictly 
positive also if we take \ai\ = \a 2 \ or d\ ■ a 2 = 1 . In particular this means that the estimate (13.41) 
is still meaningful if the angle between a\ and a 2 is 0(e p ) with < p < 1, while it clearly fails 
for p > 1. In the latter case the second atom lies inside a small cone of aperture 0(e), apex in 
the position a\ of the first atom and axis on the line joining the origin and a\. This situation 
is considered in theorem 2 where we compute the leading term of the asymptotic expansion for 
£ — > of the ionization probability. 



Theorem 2. Let us fix t > t 2 , \ai\ < \a 2 \ 
let us assume (12121) . (1233]) . (12341 . V G 5(1 



, Si • a 2 = cosxe, where Xe G [0,Xos\ 
I 3 ). Then, at the leading order for e 



Xo > 0, and 
-> 0, we have 



V(t) ~ e 2 {e- l \ ) A (e- 1 -^) Af 2 [ dxdydz 



dr} x dr} 2 F{r} x ,r} 2 \x,y,z) 
where the function F is independent of e and will be specified during the proof. 



(3.5) 



In the following remarks we briefly comment on the above results. 

Remark 3.1. The estimate ( 13.41) is valid for t larger but of the same order of r 2 , while it 
loses its meaning for t — ► oo. This is only due to the method we use for the proof, based on 
second order perturbation theory. A non-perturbative approach or a more detailed perturbative 
analysis should provide an estimate which is uniform in time. We also remark that in theorem 
2 we limit ourselves to the computation of the leading term, without making any attempt to 
estimate the remainder. Such leading term is small for e — > being proportional, as expected, 
to the solid angle that the atoms subtend at the origin. It would be interesting to extend 
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the results given here at any order in perturbation theory and in particular to verify that the 
leading term has the same behavior for e — > at all orders. 

Remark 3.2. The results in theorems 1, 2 can be understood on the basis of the original 
physical argument given by Mott, which can be described as follows. At time zero the spherical 
wave starts to propagate in the chamber and at time T\ it interacts with the atom in a\. If, 
as result of the interaction, such atom is ionized then a localized wave packet emerges from 
cii with momentum along the direction Oa\. In order to obtain also ionization of the atom in 
a 2 the localized wave packet must hit the atom in a 2 (at time r 2 ) and this can happen only if 
a 2 approximately lies on the line Oa\. It should be stressed that such physical behavior is far 
from being universal and it strongly depends on our assumptions on the physical parameters 
of the model. 



Remark 3.3. Finally we observe that our result states that one can only observe straight 
tracks in a cloud chamber. With this we do not mean that there is any focusing of the support 
of the wave packet of the a-particle along a classical straight trajectory, corresponding to the 
observed track. In fact the solution of the Schrodinger equation with Hamiltonian (12.11) and 
initial datum (12. 8p has the form 



n, r 2 , t) = F 00 (R, t)Ci(ri)C 2 (r 2 ) + J dh T c0 (R,k u t)Mn,h) ■ Ur 2 ) 
+ [dk 2 J 7 oc(R,k2,t)(i>2(r2,k2) ■ (i(n) +Jdk 1 dk 2 J r cc (R,k 1 ,k 2 ,t)cp 1 (r 1 ,k 1 )(f) 2 (r 2 ,k 2 ) (3.6) 

where the four probability amplitudes ^ooj^cOj-^Ocj^cc are localized in different regions of 
the configuration space of the whole system and therefore describe not interfering "quantum 
histories" . If one interprets double ionization as the only case of macroscopic ionization, giving 
then rise to an observable track, one can say, in a pictorial language, that is along the track 
the expected value of the position of the a particle in the states in which the track (as an 
observable) has expectation one. 

Let us outline the strategy of the proof of theorems 1, 2. The starting point is the following, 
more convenient, representation formula for the ionization probability 

A 4 t 4 N" 2 f 

V{t) = g2 dxdy 1 dy 2 \g £ 12 (x,y u y 2 ,t) + g^(x,yi,y 2 ,t)\ 2 (3.7) 
where for I, j — 1, 2, j ^ I one has 

Gfc, 2/ 2 , t) =JduJdajh j dvdt Gy (a, (3, rj, £; x, y h y 2 , t) e ^(^v,^,yuyi,t) (3 _ g) 
and dropping the parametric dependence on x,yx,y 2 ,t 
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Qij(u, a, /3, = u ■ (x + a(ar] + /3£)) — bjiij ■ n — bjd; ■ £ + CjOe + q/3 



(3.9) 



Gij(a, 0, 77, f ) = g(rj, yi)f(x + a(ai] + (3£))e 



i4>(a,f3,r),€) 



(3.10) 



<f>{a, (3,V,0 = v(v + + l {arf + fc 2 + 2an- £) 



(3.11) 



The proof of (13.71) is a long but straightforward computation and it is postponed to the appendix 
(section 7). 

Due to formula (13. 7p . we are reduced to the analysis of the two oscillatory integrals Qfj cor- 
responding to the possible graphs in the second order perturbative expansion. In particular 
Q\ 2 describes the graph in which the atom in a\ is ionized before the atom in a 2 and Q 2l the 
opposite case. Since we always assume |ai| < |a 2 |, we expect that the contribution of Q 21 is 
negligible. In fact, in section 4 we shall see that the phase 02i has no critical points and then, 
by standard integration by parts, we shall prove that the contribution of the oscillatory integral 
g e 21 is 0(e k ), for any k G N. 

The estimate of the term Q\ 2 is more delicate and we have to distinguish the non aligned and 
the aligned case. It turns out that the phase B12 has no critical points in the first case and then 
the contribution of Qf 2 is 0(e k ), for any k G N. This will be proved in section 5, concluding 
also the proof of theorem 1. 

In section 6 we consider the aligned case, where the phase B12 has a manifold of critical points 
parametrized by a vector in IR 2 . By a careful application of the stationary phase method to Q{ 2 , 
we compute the leading term of the asymptotic expansion for e — > and then we also conclude 
the proof of theorem 2. 



In this section we shall prove that the the contribution of the oscillatory integral Q 21 is negligible 
for any orientation of the unit vectors a-y, a 2 . 

Proposition 4.1. Let us fix t > r 2 , \ai\ < \a 2 \ and let us assume (12.121) . (12.131) . (12.14[) . 

V G 5(IR 3 ). Then for any fceN there exists > such that 



4. Estimate of Q 2 



21 




(4.1) 



Proof. The crucial point is that the gradient of the phase 
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©2i = u ■ (x + a(ar] + /3£)) - bjai • r\ - b 2 d 2 ■ £ + da + c 2( 5 
doesn't vanish in the integration region. To see this it is sufficient to compute 



(4.2) 



Y^p© 2 i x2 



ae 



21 



96 



facru — biai) 2 + (a(3u — b 2 a 2 ) 2 



> (act - bi) 2 + (a(5 - b 5 



a 



T-2 



(4.3) 



In the region {(a,f3) 6 M 2 | < « < 1, < (3 < a} the r.h.s. of (14.41) takes its minimum in 
(a ,/3 ), with a = j3 



2t ' 



then 



fc=i 



9% 



21 



where 



A 



h 



21 



k t 2 - n 



h \a 2 



> A 2 

— ^21 



1 - 



| ai| 



(4.4) 



(4.5) 



V2Mja V2P j 7 V l a 2| 

Notice that, under the assumptions (12.121) . (12.131) . (I2.14p and |ai| < |a 2 |, A 21 remains strictly 
larger than zero for any e > 0. The estimate ( 14.41) allows to control Q 2 \ using standard non 
stationary phase methods ([F], [Ho] . [BHJ). In fact, recalling the identity 



a e 



it 



—i div \ e 



. ib 



Vb 



a) +ie lb div 



Vb 



(4.6) 



|V6| 2 J V|V6| 2 
and performing k integration by parts we have 

Jdrjd£G 21 e^ 021 = (ie) k Jdrid£(L k G 21 )e^ @21 (4.7) 
where the operator L acts on the variables ( = (Ci> • • • j Ce) = (Vii V21 Vsi £i> £2, £3) as follows 



LG 2 \ = u 

3=1 



dG 



21 



1 86 



1 dQ ' 



21 



V ? # 21 | 2 dQ 



(4- 



and moreover 



L k G 21 = u h ..M h Dl h .. Xj G 21 (4.9) 

j'i— i*!=i 

In (14. 9 p we have denoted by D^. f . the derivative of order A; with respect to Q 1 ...Q k . From 
(TOD . (Q)l . ItOj) . (1431) we easily get the estimate 
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Ml pol r F k /■! pa r 

da dp drjd£ \L k G 21 \ < 47r— r da d(5 drjd^ 
Jo J A 21 J J J 



k rl ra 



(4.10) 

If we square ( I4.10p . integrate w.r.t. the variables x,y 1 ,y 2 and use Schwartz inequality we find 



f e 2k 

/ dxdy 1 dy 2 \<3 £ 2l \ 2 < 4vr 2 sup 

J Aoi a ,B 



d?]d^ I / dxdy\dy 2 



E ^-c^- 



ii—ife=i 



1/2' 



"I 2 



(4.11) 



From the definition of G 2 i (see (13.101) ). we have 



(4.12) 



Jl-ifc=l *1 = 1 *2 = 1 «3 = 1 



•£(|x| + oM + a|^|r 

«4=1 

The last term in (14. 12f) can be easily estimated as follows 



]T(|x| + a|7y| + a|eir 4 < £(|x + a(arj + /?£)! + 2a|r/| + 2a|C|) 

44 = 1 «4 = 1 

< £ (\/24a 2 <x + o(ot7 + /?£)>< 77 ><£>J 4 

«4 = 1 



(4.13) 



where <x> 2 = 1 + x 2 , x G M 3 . Hence 



u=i 



■ + a(m/ + £0 > k \ D *f( x + < ar l + 00)1 



12=1 



(4.14) 



«3 = 1 



Moreover, recalling the definition ( 12.261) . we get 
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6 k k 

£ \D k (ji __ 4j G 21 \ <C k <r,> k Y, \tfV{n)\ <e> fe E \D^V{i)\ 
h—jk=i »i=i »2=i 

. < x+o(m? + ^ ^ |z^ /(2;+a(m7 + \D«h(7i, yi )\ \ D i s Kc, y 2 )i 



i 3 =l 



14=1 



»5=1 



(4.15) 



Using (I4.15P in estimate (14.111) we find 



[dxdy 1 dy 2 \gi 1 \ 2 <^:cJ [drj < V > k J2\D^V( V )\ /df 1^(01 
7 I 7 ii=i 7 i 2 =i 

/to (Ej^(^)ij /^ ^X)|it/»(e,«a)i 

sup fdy [J2\K h (v,y)\ 

1 J \m=l 

J]sup ( fdy\D™h(ri,y)\ 

m=l ^ 



£ 2fc ~ 



'21 



k ||" ll^l.fc 



- 7T2fe Ck W V \\w}' k 

I—\ o 1 re 



(4.16) 



It remains to show that the last term in the r.h.s. of (14.161) is finite. From the definition (12.251) 
we have 



(2tt) 3 / 2 



dx e~ inx x^x^x™ 3 Co(x) (f>°(x,y) 



(4.17) 



where m 1 + m 2 +ra3 = m. We recall that the integral kernel (f>°(x, y) defines a bounded operator 
in L 2 (IR 3 ), with norm less or equal to one. Hence 



(4.18) 



Jdy\D™h( V ,y)\ 2 < jdx \x^xTxT Co(^)| 2 < 00 

Taking into account inequality (I4.18P in (I4.16p . we conclude the proof. 



□ 
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5. Estimate of Q\ 2 in the case d\ ■ d 2 < 1 



Following the same line of the previous section, we prove that the contribution of Q\ 2 to the 
ionization probability is negligible provided that the two unit vectors d\ and d 2 are not parallel. 
Of course, the estimate shall crucially depend on the angle between the two unit vectors. 

Proposition 5.1. Let us fix t > r 2 , a\ • a 2 < 1 and let us assume (12.121) . (12.131) . (I2.14p . 

V G 5(M 3 ). Then for any k G N there exists Ck > such that 

Jdxd yi aly 2 \gUx, yi ,y 2 )\ 2 < C k \\V\\^ k (l - ai-a^e 2 * (5.1) 

Proof. As in the case of proposition 4.1, we consider the phase 

612 = u ■ (x + a(an + /?£)) - b 2 d 2 ■ n - bihi ■ £ + c 2 a + Ci/3 (5.2) 
and we show that its gradient is strictly different from zero in the integration region. In fact 

£ (j^j 2 + (^) 2 = (aau - b 2 d 2 ) 2 + (a(3u - b^) 2 (5.3) 

The r.h.s. of (15.31) . considered as a function of the variables (a, (3), takes its minimum in 
(ai, Pt) = u ■ a 2 , ^ u ■ a x ). Hence 



^ + (j^J - ^liu- d 2 u- d 2 ) 2 + b\(ii- dxu- dx) 2 



> min{b 2 , b 2 } 2 - (u ■ d 2 ) 2 - (u ■ a^ 2 (5.4) 

Let us fix a frame of reference such that di = (0,0,1), d 2 = (sinx, 0, cosx), with x £ (0,7r]. 
Then the r.h.s. of (15. 4p . considered as a function of u = (sm8 cos (p, sin sin 0, cos (9), G [0,7r], 
G [0, 27r), takes its minimum in -u = (sin |, 0, cos |). Hence 

<9e 12 \ 2 /ae x2 



where 



fc=i 



a /o • fu u 1 A 2* h mm { «iU°2 } r, ; — fKft . 

A12 = V2 mmjbi, b 2 }Wl - cos 2 - = — V 1 - a i • a 2 (5.6) 

V 2 -H)7 7 

We notice that, under the assumptions fl 2 . 1 2 j) . (12 . 1 3j) . (12. 14j) and di • a 2 < 1, <5 i2 remains strictly 

larger than zero for any e > 0. 
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From now on the proof proceeds exactly in the same way as in the previous proposition 4.1 and 
we omit the details. 

□ 

Proof of theorem 1. From (13. 7p and taking into account propositions 4.1, 5.1 we immediately 
get the proof. 

□ 



6. THE STATIONARY CASE 



Here we consider the case ax- a 2 = 1 — 0(e q ), with q > 2. We shall see that the phase of the 
oscillatory integral Q\ 2 has stationary points when exactly this case occurs. This implies that 
the ionization probability Vit) is not negligible for e small as in the previous situation and the 
leading term of its asymptotic expansion in powers of e can be computed. 
Throughout this section we shall fix the unit vectors Si, d 2 as follows 



Si = (0,0,1), a 2 = (sinx £ ,0,cosx e ) (6.1) 

where \ e G [0, Xo £ ]j Xo > 0. Moreover, in order to characterize the asymptotic behavior of 
Qf 2 , we introduce a convenient decomposition of the unit sphere S 2 . More precisely we define 
F§ as the portion of S 2 inside a cone with apex in the origin, axis parallel to ax, aperture 8, 
< 9 < |, and we denote Tg = S 2 \ T§. For any choice of hx, a 2 as in (16.11) the corresponding 
decomposition of Q{ 2 reads 



where 



Q\ 2 = G £ 12 n + G[ 2 S (6.2) 
Q$ = fdu [da [dp [drjd£G £ 12 e? e (6.3) 

JVg JO JO J 

g[ 2 s = [du [da [dp fdrjdt G\ 2 e^ e (6.4) 
Jr s Jo Jo J 



G\ 2 = Gx2 e iS * (6.5) 

sinx e 1 - COSXe u fa R ^ 

Oe = — b 2 % H b 2 r] 3 (6.6) 

= u ■ (x + a(ai] + /3C,)) - bi^ 3 - b 2 r] 3 + c 2 a + CxP (6.7) 
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We shall analyze the asymptotic behavior of the two oscillatory integrals Q £ { 2 and Ql 2 separately. 
We first show that the phase has no stationary points in T§ and then the contribution of Q^ 2 
is negligible. 

Proposition 6.1. Let us fix t > r%, a\, a 2 as in (16. ip and let us assume (" 12. 12j) . fl 2 . 1 3 j) . (I2.14p . 

V G 5(1R 3 ). Then for any k G N there exists > such that 

Jdxd yi dy 2 \gi> 2 n (x,y u y 2 )\ 2 < C k (-^) 2 " (6-8) 

Proof. If we denote u = (sin 9 cos 0, sin 9 sin (ft, cos 9) G T§, we have 




= a 2 (a 2 + f3 2 ) + b 2 + b 2 , - 2a(bi/3 + b 2 a) cosfl 



> a 2 (a 2 + /3 2 ) + b 2 + b 2 - 2a(bi/9 + b 2 a) cos (6.9) 
The r.h.s. of (16. 9p takes its minimum in (a 2 , /3 2 ) = ( „ cos#, ^ cos#), then 




(6.10) 



where 

A = V2 minfox, b 2 }sinfl = A min {M, H> sin g (6 .n) 

P 7 7 

Exploiting estimate (16.101) . (16.111) it is now straightforward to obtain (16. 8p . 

□ 

Remark 6.1. We notice that the estimate (I6.8P is still meaningful if we choose the angle 9 
depending on e and such that 9 = 0(e d ), with < d < 1. This in particular means that only a 
small fraction of the unit sphere, of area 0(e 2 ) around the direction Si, can give a non trivial 
contribution to the ionization probability. 

Let us consider the oscillatory integral Ql 2 . ft turns out that the phase has a manifold of 
critical points in the integration region, parametrized by a vector in R 2 . Therefore we fix the 
variables (?7i,^ 2 ) G M 2 as parameters and we write Q\ 2 i n the form 



l £ (Vi,V2)= dqG £ 12 (q]ri 1 ,n 2 )e 



(6.12) 
(6.13) 
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where 

fl={q = («, a, /?,7/3,£) \ueTg, ae [0,1], /? G [0,a], r/ 3 G R, ^Gl 3 } (6.14) 

In the next lemma we show that for each value of the parameters (771,772) the phase in (16.131) 
has one, non degenerate stationary point. It is relevant that the value of the phase and of the 
hessian of the phase at the critical point do not depend on (771, rf 2 ). 

Lemma 6.2. For each (r/i, 7/2) G ^ 2 the phase Q(q; 771, 772), q G Q, has exactly one critical point 



where 



q = (u", a", (3\ rf z , & & £3°) (6-15) 



«° = (0,0,1), a° = ^, = 7/° = -^, (6.16) 

a a a 

o_ si + Mi _ ^ + M2 _ Ci 

Cl — 7 ) S2 — 7 ! S3 — ~ I - 1 ' J 

bi bi 



and moreover 



6° = B(q ;Vi,V2) = x 3 + 



bid , b 2 c 2 



|L> 2 e°i = |D 2 e(go;r/i,r/ 2 )| = a 4 bt 



(6.18) 
(6.19) 



Proof. In order to compute the critical points of the phase (16.71) as a function of q G Q it 
is convenient to write u = (fx, v, \f\ — fx 2 — 1/ 2 ), where (fx,v) G IR 2 with fx 2 + v 2 < sin 2 
Therefore 



©(<?; = fxw! + vw 2 + a/1 - fx 2 - v 2 w 3 - b 2 77 3 - bi^ 3 + c 2 a + Ci/3 (6.20) 

where we have introduced the short hand notation 

«; = («Jx, w 2 , w 3 ), wj = xj + a(arfj + (3£j) (6.21) 



By an explicit computation, one finds that the critical points are solutions of the system 
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•96 fj,w 3 

— — = wi = = 

dli y/l-IJ?- v 2 

(96 vw 3 

w 2 



dO 



, - aurji + aisr)2 + a\/l - fi 2 - u 2 r] 3 + c 2 = 

— = + 01/6 + Oa/1 - yLi 2 - + Ci = 

56 



ay 1 — jj, 2 — v 2 a — b 2 

= 
oi//3 = 



— = o v / l-/i 2 -^,3-b 1 = 



First we notice that a and (3 cannot be zero, otherwise from (I6.26p . (16.291) one would have 
b 2 = bi = 0. Then from (IOTP . (I6T28D we have /x = v = and from fl6\26|) . f[6T2?j) we have 
a = P = ^. Exploiting the remaining equations it is now trivial to find the unique solution 
f)6.16p . (16.1 7p . Furthermore the value of the phase at the critical point (I6.18P is easily obtained. 
For the proof of (16.191) we need the second derivatives of the phase evaluated at the critical 
point 



d 2 e _ 


-w 3 , 


d 2 e 

dfidu 


= 0, 


o 2 e 

p> ^ = 


d 2 Q 




d 2 Q 


= 0, 


d 2 Q 




= bi, 




d[xd£ 3 


d 2 e 

dv 2 


-w 3 , 


d 2 e 

dvda 


= ar] 2 , 


d 2 e 

dvdp 


d 2 e 


= 0, 


d 2 e 


= bi, 


d 2 e 






dvd£ 3 



a6 



o 2 e 

dud/3 



a6 



d 2 e 

dh>dr] 3 



d 2 e 

dfidr] 3 
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Q2Q ^ Q2Q Q2 Q ^ Q2 Q 

da 2 ' dad(3 ' dadr]^ ' dad£j 

Q2Q _ Q2Q Q2Q Q2Q Q2Q 



d(3 2 ' d(3dr] 3 ' d(3d^ ' d[3d& ' d(3d& 

d 2 e d 2 e d 2 e 

The computation of the hessian is now a tedious but straightforward exercise and it is omitted 
for the sake of brevity. 

□ 



We are now ready to conclude the proof of theorem 2. 

Proof of theorem 2. Exploiting the stationary phase theorem ([F], |Hoj . |BHj ) and the 
previous lemma, we find the leading term of the asymptotic expansion of (16.131) for e — > 

I%Vi,V2) ~ e" °G 12 (qo;Vi,V2)e l5o e^ (6.31) 

where //q denotes the signature of the hessian matrix at the critical point and moreover 

Sq = — lim Sm ^ £ b 2 r/i (6.32) 

e^O £ 

In particular S ^ if Xe — 0(e) and S = if % e = 0(e 9 ), g > 1, or Xe — 0. From (I6.12p we 
also obtain 

GI2 ~ e* e § y^dTfe G 12 (g ; r/x, 77 2 ) e ia ° (6.33) 

We notice that the integrand in (16.331) is a function of x (position of the a-particle), yj (mo- 
mentum of the j'-th ionized atom) and 771 , 772 . Hence we denote 

F(rj 1 ,r ]2 ;x,y 1 ,y 2 ) = (2vr) 4 G 12 (g ; e i5 ° (6.34) 

and, taking into account (13.71) . proposition 4.1, (I6.2p . proposition 6.2 and (16.331) . we find 



drfidri2F(Tii,Ti2]x,y,z) 



2 



(6.35) 



Using the definition of a, bi, 7~i in (16.351) . we easily get formula (13.51) . It remains to show that 
the integral in (I6.35j) is finite. From the definitions (I3.10p . (12.261) . (12.251) and the boundedness 
of h(£, y) we have 

\F( Vll r, 2 ;x,y,z)\ < c\V( Vl ,r ]2l r ) 3 )\\V(e i ,e2,e3)\\f(wo)\ (6.36) 
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where 



M 



m Oi Ti ai Ti 



r 2 



-i 7 



r 2 



M 2 M l fl ll / \ 2 M I«2| . . 

4 3 = -e—u{y) , w = x 3 - e u{y) - e u{z) 

m m 7 ■ ■ 



m 7 



and u)(y) is defined in f 1 2 . 2 3 j) . Then we write 



dxdydz 



dr) 1 d7) 2 F(7)i,7]2;x,y,z) 



Kcjdx^dydzlJdvMVi^^^Wie^^^lj jdx 3 \f(w )\ 2 
Using the Schwartz inequality in the integral with respect to (771 , 772) we have 



(6.37) 



(6.38) 



dxdydz J d7)idr]2 F(t)i, 7)2,x, y, z) 

<cjdzjd Vl d V2 <7 ?1 > 4 <r / 2> 4 |y(r ?1 ,r /2 ,r /3 )| 2 J dy J dxrdx^V^^il)? J dx 3 \}\w )\ 2 

(6.39) 

The last integral is finite due to the assumptions on V and this concludes the proof of the 
theorem. 

□ 



7. APPENDIX 



Here we give a proof of the representation formula (13.71) . The relevant object to compute is the 
probability amplitude in (13.31) 



T{R,kx,k 2 ,t) = dndr 2 (f}i(r 1 ,k 1 )(f)2{r2,k 2 )if2{R,ri,r 2 ,t) 



We notice that 



ei tH »H ie -i tH ° = Wx(t) + W 2 {t) 
Wj(t) = e^e^Vje-^e-^ 
where Vj denotes the multiplication operator by 



(7.1) 



(7.2) 
(7.3) 
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V 3 {R,r ) = V{T\R-r J )) (7.4) 

Let us rewrite the r.h.s. of (17. ip in a more convenient form. We observe that the operator 
Wj{t) acts non trivially only on the variable R and rj. Exploiting this fact it is easily seen that 

T = Fi2 + Fix (7.5) 
F lj {RMMt) = —jdt l jdt 2 Jdr l dr 2 ^ {W^W^m^R.r^) (7.6) 

where l,j = 1,2, j 7^ I. Furthermore, using the specific factorized form of the initial state of 
the system, we have 

e^V^-M^^yR) (7.7) 
where E(k) = and Vj(-, kj), j = 1, 2, is the multiplication operator by 



Vj(R, %) = dr </>Jr, kjW-y-^R - r))Q(r) (7i 



The r.h.s. of (17.81) can be more conveniently written in terms of the Fourier transform V of the 
interaction potential as follows 

Vj(R, kj) = ^ 2 Jdi J^^X g(£, 7%) (7.9) 

where g(£, y) has been defined in (I2.26p . Exploiting (17. 9p and the explicit expression of the free 
propagator we have 



and 



^' 2 ^°^(-,A; i )e-^ 2X V)( J R) 

j 3 ^g(tlk) e^ e+ ^~^U(R + ^) (7.10) 



= ' ' kr "> fjd^gi^k^^h) 

. e i^ 2+ T^^< + ^<) e<-^- a i<)^(R+^ v+ ^f) (7 11) 

r V M7 ' M7V v 
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Finally we consider the time-dependent phase factor in (17.71) . We notice that 

|(£(%) + \E \) + j(E(k) + \E \) = ^ (hwfrkj) + t 2 w(lh)) (7.12) 

Taking into account (17.111) . (I7.12p . (12.91) . and rescaling the time variables according to t\ = ta, 
t 2 = t/3, we can rewrite (17. 7p as follows 

Tij{R, h, k 2 ,t) 



ft 2 e y<?2 j j j v 7 M 7 2 

(7.13) 

We observe that 

^ = -0(1) (7.14) 

and this means that the first exponential in the integral in (I7.13P has a slowly oscillating phase 
for e <C 1. On the other hand 

'" M a = 0( £ -) (7.15) 



r«7 2 m 

and therefore the last exponential in the integral in (I7.13P has a rapidly oscillating phase for 
e <C 1. Denoting R = jx, kj = 7 _1 ?/j and using the notation (12.241) . we find 

/ ' j^ e - ik ^- ik ^ du da d/3 dCd V g( V , yj )g(C,yi)f(x + a(ar, + M)) 



h 2 e Js 2 Jo Jo 

e i[z-(»7+£)+f (c«, 2 +/3£ 2 +2ar,-£)] g f [«-x+«-a(a ?? +/30-b J <i J -r)-b ; a r C+c J a+c i /3] 

s ^^e-^-riV^^,^^^) (7.16) 

where in the last line we have used fl3~9D . fl3~T0|) . flBTTTj) . From (Q, (JUJ) and (1771731) we 
obtain 

V(t) = 7 -3 jdxdy x dy 2 (^(t^ 7~V, 7^2/2) + T 21 {^x, j~ x yi, T 1 y2) | 2 

Jdxd yi dy 2 \gi 2 (x, yi ,y 2 ) + ^(x, yi, ?/ 2 )| 2 (7.17) 



A 4 t 4 N" 1 



h A e 2 

and this concludes the proof of (13.71) . 
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